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ABSTRACT 

We consider the effect of planetary spin on the planetary radial velocity (PRV) in dayside spectra of 
exoplanets. To understand the spin effect qualitatively, we derive an analytic formula of the intensity- 
weighted radial velocity from planetary surface on the following assumptions: 1) constant and solid 
rotation without precession, 2) stable and uniform distribution of molecules/atoms , 3) emission 
models from dayside hemisphere, and 4) a circular orbit. On these assumptions, we find that the 
curve of the PRV is distorted by the planetary spin and this anomaly is characterized by spin radial 
velocity at equator and a projected angle on a celestial plane between the spin axis and the axis of 
orbital motion X p in a manner analogous to the Rossiter-McLaughlin effect. The latter can constrain 
the planetary obliquity. Creating mock PRV data with 3 km/s accuracy, we demonstrate how A p and 
the spin radial velocity at equator are estimated. We find that the stringent constraint of eccentricity 
is crucial to detect the spin effect. Though our formula is still qualitative, we conclude that the PRV 
in the dayside spectra will be a powerful means for constraining the planetary spin. 

Subject headings: planets and satellites: fundamental parameters - techniques: radial velocities - 
techniques: spectroscopic 



1. INTRODUCTION 



Planetary spin is one of cruci al factors that gov- 
ern t h e climate of exoplan et s (e.g . IWilliams fc Kastind 
19971 IWilliams fe Pollardl [2001 IHeller et al. "Una 
Cowan et al.ll2012t ) and has a potential to constrain plan- 



etary formation theor y (e.g. lAgnor et alll999l iChambersI 
120011 : iKokubo fe fdal I2007t ). Photometric variation of 
a planet will enab le us to estimate t he rotation period 
(iPalle et ap 120081) and the obliquity (|Kawahara fc Fuiiil 
120101 1201H iFuiii fc Kawaharall2012l ) in the near future. 
However these methods are applicable only for a planet 
having significant inhomogeneous surface, such as coex- 
istence of ocean and lands or clouds, and also require a 
long-term observation with a sophisticated instrument of 
direct imaging. 

Recently planetary radial velocity (PRV) of the non- 
transiting planet, r Bootis b has been measured with 
the aid of the carb on monoxide absorption in the therma l 
dayside spectrum (jBrogi et al.ll2012t iRodler et al~ll2012[ ) . 
iBrogi et "all <|2012f ) detected the change in the radial com- 
ponent of the planet's orbital motion and obtained the 
semiamplitude of the PRV, K p = 110.0 ± 3.2 km/s. 
Combining the stellar radial velocimetry with it, they 
evaluated the orbital inclination and mass of r Bootis 
b. Though the PRV curve is primarily dictated by the 
planet's orbital motion, the planetary spin also has pos- 
sible effect on the PRV in the dayside spectrum. Gas 
giants in the solar system have considerable spin velocity 
at equator, 12 km/s for Jupiter and 10 km/s for Saturn. 
Even for most hot Jupiters which are likely to be in a 
synchronous orbit, are expected to have non negligible 
spin velocity driven by their rapid orbital motion. For 
instance, WASP 19b will have ~ 9 km/s of the spin ve- 
locity at equator if it is tidally locked. The aim of the 
paper is to develop the method to derive the planetary 



spin-orbit alignment and the spin velocity from time se- 
ries analysis of the PRV. This concept can be explained 
by an analogy to the R o ssiter-McLaughlin effect (RM 
effect lOueloz et al.l 120001 IQhta et all 120051 IWinn et al.l 
l2lMlGaudi fc Winnll2007 . and references therein). The 
RM effect is an anomaly of Stellar radial velocity (SRV) 
caused by sequent occultation of a rotating stellar disk 
by a transiting planet, and is used to measure the pro- 
j ected angle of th e orbital axis and the stellar spin axis 
(jOhta et al.ll2005l ). Likewise, non-uniform emission from 
a planet, which is generally stronger near the sub-stellar 
direction, induces an anomaly in a time series of the PRV. 
In this paper, we demonstrate how the planetary spin af- 
fects the PRV assuming a simple solid rotation of a planet 
and derive an analytic formula of the PRV anomaly with 
simple intensity distribution models. 

2. METHODS 

We divide the PRV into the radial velocity components 
of the planetary center system iv,orb(@) and the plane- 
tary spin U r ,rot(Q)i 



V r {&) =lV,orb(0) + tV,rot(0)- 



(1) 



On the assumption of a circular orbit for simplicity, the 
PRV by the orbital motion is expressed as 



H r ,orb(0)=-KpCOs6, 



(2) 
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where K p is semiamplitude of the radial velocity by the 
orbital motion and O is the orbital phase. We define 
the phase angle between the line of sight and star-planet 
direction a, 

cos a = eg • eo = sin i sin O 

( tt/2 - i < a < tt/2 + i ), (3) 

where e$ = (cos a, sin a, 0) T and eo = (1,0, 0) T are unit 
vectors from the planetary center to the star and the 
observer and i is the orbital inclination (Figure [T] A) . 
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Figure 1. Geometric configuration. Panel A: the spherical coordinates on a planet. The stellar direction is embedded in the x-y plane. 
The unit vectors eo, es, and, en are explained in the text. Panel B and C presents schematic pictures from observer's point of view 
(y-z plane). In panel B, the blue shifted emission dominates on the dayside hemisphere (drawn by white), which causes the blue-shifted 
anomaly in the PRV. The planetary spin orbit misalignment X p is defined as the projected angle between the spin axis and the orbital axis 
on the celestial plane (Panel C). Panel D indicates relation of the angles, the obliquity £, the orbital inclination i, the spin inclination i p , 
and X p . The uncertainty range of i p du e to the degeneracy between sini p and ui apln R p for a solid rotator is shown. Panel E: possible range 
of obliquity ( as a function of A p (Eq. [10]). Each shade corresponds to dilferent orbital inclination i. 



Since the observed absorption consists of ensemble of a 
Doppler-shifted line from each position on planetary sur- 
face, iv,rot(@) depends on the distribution of absorption 
lines on the surface. We use the intensity distribution 
instead of the distribution of absorpt ion lines assuming 
the uniform distribution of molecules. lOhta et ail (|2005l ) 
showed that the intensity-weighted velocity is in agree- 
ment with the Doppler shift, v riT0 t/c = Ais/v on the as- 
sumption that the frequency shift is much smaller than 
the line frequency (Eq. [20] in their paper). Though 
the precise value measured by real observation will de- 
pends on details of methods and instruments, the aim 
of the paper is to qualitatively understand the behav- 
ior of iv,rot(©)- Hence we regard the intensity- weighted 
velocity, 

v rjot = -L- f dnW(ct>,9;a)V r (<j),6) = (V r ), (4) 



/ dnW(<j>,6;a) 



(5) 

as the PRV of the spin, iv,rot(9), where the ^(a) is the 
phase function (total intensity of emission), W(4>,6;a) 
is the intensity distribution normalized so that ^(0) 
is unity, V r ((f>, 6) is the radial velocity of a facet on 
surface measured on the planetary center system, and 
dfl = sin 6d0d<j) (sec Fig. [TJV). 

A naive expectation is that W(4>, 9; a) has higher value 
at the dayside hemisphere than that at the nightsidc 
hemisphere (Fig[T]B and C). This is absolutely the case 
for the scattered light. It is also likely to be the case 



for thermal emission on the assumptions that the energy 
injection is mainly attributed to incident light from the 
host star. 

In this paper, we regard a planet as a solid rotator and 
no precession, 

V r (4>,0;T) = K I0 t(— sin r sin <j> sin 6 + cost cos 6), (6) 

where t = Q + ir/2 — X p is the projected rotation angle 
between the spin axis and the stellar direction on ce- 
lestial plane and the projected orbital phase angle 9 is 
expressed as 



tan 9 = cosi tan 9. 



(7) 



The maximum radial velocity of the solid rotator K TOt 
represents the strength of the spin, 



K Iot = 27rw spin i?p sinip, 



(8) 



where i p is the spin inclination. Thus the degeneracy 
between smi p and u; sp i n is inevitable for a solid rotator. 
The planetary obliquity £ is expressed as, 



cosC 



cos i cos % r , 



sini sinip cos X p , 



(9) 



where we introduce the planetary spin orbit misalign- 
ment Ap, defined by the projected angle between the spin 
axis and the orbital axis on the celestial plane. Figure Q] 
D displays a schematic picture of these angles. Though 
i p is difficult to know directly, if knowing X p from the 
PRV time series analysis, one can constrain the obliq- 
uity £ for a given i. Figure [T]E shows the possible region 
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Figure 2. Panel a: cos© as a function of ©. Panels b and c: 
F(0) (black) and phase function >1>(0) (red) for the IRR and DR 
models. Two vertical lines correspond to SC and IC. 

of £ for given A p , 

sin^ 1 (sinz| sin A p |) < £ < n — i 
for < A p < tt/2, 37r/2 < X p < 2tt, 



i < C < n ~ s i n 1 (sin i I sin I ) 



(10) 



for tt/2 < A p < 3tt/2. 

If the intensity distribution is symmetric about the 
stellar direction (W(<j),9;a) = W{<j), tt/2 - 9; a)), the 
antisymmetric cos 9 term in equation ^ vanishes and 
equation ^ reduces to 



*V,rot = -K Iot F(a) cos (6 + Ap), 



(11) 



where 



F(a) 



*(a) 



dQ.smesm(t)W(<f>,6;a). (12) 



The </>-symmetric component of W{<$>, 9; a) does not con- 
tribute to tv.rot. In equation (|11[) . one can interpret that 
F(a) represents the inhqmogeneous pattern of the plan- 
etary surface. The cos (O + A p ) term is due to apparent 
change of daysidc along orbit and does not depend on 
the intensity distribution. As increasing i, the deviation 
by this term from the cosine curve of the orbital motion 
becomes stronger as shown in Figure [5] a. 

3. ANALYTIC FORMULA OF THE SPIN EFFECT WITH 
RADIATION MODELS 

Considering two specific models of the intensity 
distribution, the instant re-radiation (IRR) model 
and the dayside re-di s tribut io n (DR) mode l (e.g . 
Lopez-Morales fc Seaeen 120071 : iCowan fc Ago] 120111 : 
Smith et al.l 120121 ). we derive an analytic expression of 
the spin effect. 

The IRR model assumes the input energy from a host 
star is instantly and isotopically emitted from the plan- 
etary surface. The isotropic reflection for the scattered 



light or no redistribution of energy around surface for the 
thermal emission satisfies this assumption. The intensity 
distribution of the IRR model is expressed as, 



— (e s • e R )(e R ■ e ) 

27T 

for -tt/2 + a < <f> < tt/2, 
for elsewhere, 



(13) 



and the phase function 

^iRR( a ) = —(sin a + (n — a) cos a), 



where e R = (cos <fi sin 9, sin 4> sin 9, cos ( 
(fT2l. we obtain 



(14) 

Using equation 



FmR(a) 



37Tsina(cos a + 1) 
16[(7T — a) cos a + sin a] 



(15) 



The DR model assumes rapid energy redistribution on 
the daysidc. 



[ for elsei, 

* D R.(a)=cos 2 (|) , 

fl»( )-±^(H) 



a < (j) < tt/^ 16 ^ 
ir elsewhere, 

(17) 
(18) 

Figure [2}3 and c display F and \1> for both models as 
a function of 9. Substituting FmR(a) (F^(a)) into 
equation (jllj). we obtain the analytic expression of w rjro t- 
Figure [3] (left) shows zv jro t for different A p and i. One 
can see characteristic features of the PRV depending on 
A p and i. By fitting this anomaly, we can estimate X p 
and K Iot as will be demonstrated in next section. The 
difference of v r . ro t between these models is not significant. 

For the exact edge-on orbit, equation (fTT|) reduces to 



v r ,rot = ±27rw S pi n i?pCosC-F(a), 



(19) 



where a positive sign is for 7r/2 < < 37r/2 and nega- 
tive for 6 < 7r/2 or > 37r/2. Equation (|T9| indicates 
that the phase information of the term cos (6 + A p ) dis- 
appears except for its sign in the edge on limit. However 
one can know whether the spin is prograde or retrograde 
against the orbit. 

While we have considered the line shift so far, the spin 
rotation also induces the line broadening. We simply 
estimate the doppler broadening by considering the vari- 
ance, 



(V?) (V r f 



(20) 



where the definition of () is same as that in equation @ . 
The right panels in Figure |3] show for i = 45°, 85° 
and 0°. Typical broadening width is ~ 0.4 K rot . The 
0-dependence is smaller than that of the line shift. The 
dependence of the emission models on <tl is not signif- 
icant. If using extremely high dispersion spectroscopy 
(for instance, R = 300,000 for K Iot = 10 km/s), one 
may derive this characteristic feature of the broadening 
variation. 

4. DEMONSTRATION 

We demonstrate the spin effect by mocking the PRV 
time series. Figure [4] shows the mock PRV curves with 
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Figure 3. Left panels show the PRV anomaly by planetary spin as a function of © in an inclined orbit (i = 45°; top) and a typical 
transiting system (i = 85°; middle), and a face-on orbit (i = 0°; bottom). Each color corresponds to different X p . Solid lines assumes the 
IRR model, while dotted lines indicate the DR model. Right: line broadening due to a solid rotation for different inclination (same as the 
right panel). Standard deviation normalized by K ro t is shown by solid (dashed) lines for the IRR (DR) model. 



the IRR model and its fitting curves. We set a = 3 
km/s precision, roughly corresponding to the accuracy 
obtained when lines are barely detected and resolved 
with R ~ 100, 000, though it can be improved by increas- 
ing sensitivity and the number of lines. The top panels in 
Figure [4] assumes a non-transiting system with i = 45°, 
K p = 100 km/s, K Tot = 10 km/s, and X p = 45° (left; 
case A) or 0° (right; case B). We take 100 data points 
and avoid the ±30° range around the inferior conjunc- 
tion (IC) since the planetary signal in this range declines 
below 7 % (20 %) of the maximum value at 9 = 90° for 
the IRR (DR) model (see Figure 2bc red lines). We use a 
Gaussian prior for the stellar mass with typical 5 % un- 
certainty and simultaneously fit the stellar mass, i, X p , 
and -Krot using a Markov chain Monte Carlo (MCMC) 
algorithm. The estimated inclination is well constrained 
i = 45 ± 3° since it is almost determined by the whole 
amplitude. The spin effect can be detected above 3 a for 
these cases; K rot — 10 ± 2 km/s with X p = 40^ ° (case 
A) and K mt = 8 ± 3 km/s with A p = 3 ± 7° (case B). 

Though we have assumed an exact circular orbit so 
far, uncertainty of small eccentricity, e can be a possible 
source of systematics. For e <C 1, one can express the 
velocity modulation due to e as, 

v r ,orb ~ A" p [cos 9 -(- ecos (29 — w)], (21) 

where w is the argument of periastron. The modulation 
has the amplitude of eK p and double frequency. We per- 
form the MCMC leaving e and u> as fitting parameters 
with equations (fTT|) and (|2"Tj). Assuming a typical con- 
straint of e < 0.02, we obtain slightly worse constraints, 
K mt = 8±3 km/s and \ p = 45±f ° for the case A. 
For the case B, the eccentricity uncertainty makes the 
spin effect to be unconstrained (1 a detection level of 
A" ro t) due to the correlation between e and K mt . Thus 



e < 0.02 is marginal to detect the spin effect since the 
uncertainty of the modulation is comparable to the am- 
plitude of the spin velocity, that is, K p Ae ~ K mt . If 
assuming e < 0.002, which satisfies Ae <C K Iot /K p , we 
obtain A rot = km/s and 2+g ° even for the case B. 
As shown in the top left panel, u r ,rot for the case B is 
resemble to the eccentricity modulation (Eq. [21] ) in its 
curve. Therefore the case B is more sensitive to the e un- 
certainty than the case A. A more stringent constraint 
of e than K TOt /K p is important to detect the spin effect. 

For the time being, a realistic target of the PRV mea- 
surement will be confined to hot Jupiters. We also con- 
sider the application of the precise PRV measurement to 
hot Jupiters. Orbital period of hot Jupiters is generally 
short (P ~ a few days, typically). In a synchronous orbit, 
the rotation period equals to the orbital period. Contri- 
bution of the orbital motion to zv,rot is 2-kR p /P = 1 — 3 
km/s for typical hot jupiters and reaches 2irR p /P = 9 
km/s for the extre me case, WASP 19 b (P = 1.2 day 
and R p = 1.39E?: iHellier et alJl20ll . We create the 
mock PRV of a tidally locked transiting planet with pa- 
rameters of WASP 19b (A rot = 228 km/s and i = 79°). 
Since i is known for the transiting planet, we assume 5 
% uncertainty of the stellar mass again. We fit the PRV 
curve avoiding the range ±45° around transit, in which 
the planetary signal declines below 5 % (IRR) and 15 % 
(DR) of the maximum. We also take the eccentricity un- 
certainty into consideration by adopting e = 0.046 and 
uj = 0° for input. Since current uncertai nty of WASP 
19 b is e = 0.0046°:E] 4 4 28 and u = 3 ± 70° (IHellier et al.1 
120 111 ), we assume a Gaussian prior with a e = 0.004 and 
er w = 70°. We obtain a marginal detection of the spin, 
K Tot = ll±l km/s and X p = 4±^ . The fitting with 
more stringent constraints a e = 0.001 and a u = 20° pro- 
vides A: ro t = 13 ±3 km/s and X p = 11±^ °. The bottom 
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Figure 4. Mock data of the PRV and residuals with 3 km/s accuracy. Top panels: mock PRV curves for oblique (left; \ v = 45°) and 
aligned (right; X p = 0°) planets for a non- transiting system (i = 45°). Black and red curves in the upper subpanels are the model fit 
and cosine fit of the data. The bottom subpanels display the residual of the model fitting (black) and the cosine fitting (red). Solid 
curves indicate the model fitting minus the cosine fitting, that is, the anomaly from the cosine curve. Bottom panels: simulated data for 
a transiting system mocking the WASP 19 b-like system (i = 79°). Bottom left is a tidally locked case (the planet's rotation is directly 
driven by the orbital motion). Bottom right: an aligned planet but having the rapid retrograde spin (X p = 180°) with K Ia t = 15 km/s. 



left panel shows the mock PRV and its fitting results for 
the latter. 

Most hot Jupiters are likely to be in a synchronous 
orbit with £ = (therefore A p = 0) due to tidal force. 
However, this expectation has never been proved obscr- 
vationally. In particular, the assumption of the tidal lock 
is not obvious for a planet in a highly eccentr ic orbit such 
as HD 80606 b (e = 0.93; lFossev 'eTaI|[2009t) . Finally we 
consider an extreme planet with a rapid retrograde spin 
C = 180° and K rot = 15 km/s. We note that retrograde 
here means the retrograde rotation of the planetary spin 
( not the stellar spin ) against the orbital revolution, like 
Venus. Performing the same fitting process as the tidally 
locked case (with Ae = 0.004 and Aui = 70°), we obtain 
-Krot = 17 ±3 km/s and A p = 181 ±17° with a character- 
istic feature of the residual as shown in the bottom right 
panel. 

In this section, we excluded the light curve around the 
IC where the signal is weaker as is clear from ^(6) in 
Figure^] As shown in the bottom subpanels of Figure 0] 
most characteristic features appear around the IC though 
it depends on i (see also Figure [3]) . This is one of main 



difficulties to detect the spin effect practically. Though 
current detection of the PRV is far away from the IC 
(0.5 < 9 < 2.5 for iBrogi et al.ll2012D . future detection 
in the outer range is of importance to detect the spin 
effect. 

5. DISCUSSION AND CONCLUSIONS 

In this paper, we have estimated the observable shift 
of the planetary absorption by substituting the inten- 
sity weighted radial velocity (Eq. [4]). It is known 
that the cross-correlation method used in the RM ef- 
fect has a few - 50% systematics of the amplitude of the 
velocity anomaly when using the intensity weight (e.g . 
iWinn et al.|[2005l: iTriaud et al1[2009t iHirano et al1l2010l) . 
Moreover, even two cross-correlation methods used in dif- 
ferent instruments make < 30 % difference in the am- 
plitude of the RM effect (T. Hirano in private commu- 
nication). Hence more sophisticated modeling adapted 
to the details of measurement wi ll be needed to ob - 
tain precise estimates practically. IHirano et al.l ([2010D 
found that the radial velocity anomaly obtained by the 
cross-correlation method is larger than that predicted 
by the intensity weight method and that the bias tends 
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to be larger as increasing the spin velocity. It makes 
detection of the spin effect easier. Detailed structures 
of planet ary surface such as u neven molecular distribu- 
tion (e.g. iBurrows et al.ll2010( ) or winds may affect the 
PRV in principle. The PRV in transmission spectra 
([Snellen et al.l l2010h might help to resolve a degen eracy 
between winds and spin (see also lSpiegel et al.ll2007t ). We 
postpone these problems until a next paper. 

In conclusion, we have shown how the planet's rotation 
affects the PRV in the dayside spectra of exoplanets. We 
found that the spin effect of the solid rotation on the 
PRV is characterized by the projected angle between the 
orbital axis and the spin axis, X p . We also showed that 
the precise measurement of the PRV enable us to con- 
strain the planetary obliquity via X p and the spin period 
via If rot . 
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